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Abstract 


This report contains a description of the test facilities and software utilized during a 
joint NASA! aerospace industry study of improved control laws and desired inceptor char- 
acteristics for a candidate supersonic transport aircraft design. Details concerning the 
characteristics of the simulation cockpit, image generator and display systems, and motion 
platform are described. Depictions of the various display formats are included. The test 
schedule, session log, and flight cards describing the maneuvers performed are included. A 
brief summary of highlights of the study is given. Modifications made to the industry- 
provided simulation model are described. This report is intended to sene as a reference 
document for industry researchers. 


Introduction 

A joint NASA/industry team has been examining the 
feasibility of producing and operating a profitable supersonic 
transport aircraft under sponsorship of the High-Speed Re- 
search (HSR) Program. This report describes one such study, 
an investigation by the Guidance and Flight Controls Inte- 
grated Technology Development Team, conducted in the 
spring of 1998 at NASA Langley Research Center (LaRC). 
This study utilized the Langley Visual/Motion Simulator 
(VMS) to develop control laws to mitigate the effect of struc- 
tural modes on passenger ride and crew flying qualities. The 
study also looked at ways of quantifying the effect of vary- 
ing lateral inceptor authority on flying qualities. This study 
was designated LaRC. 4 and was the fourth experiment at 
LaRC under the HSR program. 

The intent of this report is to serve as a reference for the 
industry investigators of the separate elements of the LaRC.4 
study when reporting on specific results of those elements. 


Abbreviations 


ANSI 

American National Standards Institute 

CRT 

cathode ray tube 

DASE 

dynamic aeroservoelastic 

HSR 

High-Speed Research 

IMC 

instrument meteorological conditions 

LaRC 

Langley Research Center 

LaRC.4 

The fourth in a series of planned HSR 
tests at Langley Research Center 

SMC 

structural mode control 

SSL 

Secure Sockets Layer 

TCP/IP 

Transfer Control Protocol/Intemet 
Protocol 

VMS 

Visual/Motion Simulator 

XVS 

eXtemal Vision System 


Facility description 

Host processor 

The host processor used for the LaRC.4 experiment was 
a Convex 3840, running the vehicle model at a major frame 
time of 16 ms (62.5 Hz), with some components (actuators 
and dynamic aeroservoelastic models) running at higher 
frame rates (usually 5 times per major frame, or 3 .2 ms [3 1 2.5 
Hz]). Cockpit, motion, and visual input/output were per- 
formed at 32 ms (31.25 Hz). 

Motion base and cueing algorithm 

The motion base utilized in LaRC.4 was the Langley 
Visual/Motion Simulator (VMS), which was designed and 
manufactured by The Singer Company in 1972. It is de- 
picted in figure 1 . The VMS is a synergistic six-degree-of- 
freedom motion platform with the nominal amplitude, rate, 
and acceleration limits given in table 1 . 

The motion cueing algorithm utilized was a modifica- 
tion of normal VMS motion cues. Normal VMS motion 
cueing uses coordinated adaptive washout filters as described 
in references 1 and 2. The modification to these conven- 



Figure 1 . LaRC VMS exterior view. 
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Figure 2. Motion cueing compensation for DASE. 


Table 1 . LaRC VMS Static Characteristics 


Axis 

Position 

Limits 

Velocity 

Limits 

Acceleration 

Limits 

Surge 

49 in. forward 
48 in aft 

+24 in./sec 

+0.6 g 

Sway 

±48 in. 

+24 in./sec 

±0.6 g 

Heave 

39 in. up 
30 in. down 

+24 in./sec 

±0.8 g 

Roll 

+22 deg 

+15 deg/sec 

±50 deg/sec 2 

Pitch 

30 deg nose up 
20 deg nose down 

+ 15 deg/sec 

±50 deg/sec 2 

Yaw 

+32 deg 

+ 15 deg/sec 

±50 deg/sec 2 



Figure 3. LaRC VMS interior view. 


tional motion cues was to take that portion of cockpit accel- 
erations due to dynamic aeroservoelastic behavior and add 
them directly to the output of the washout filters to reduce 
phase lag in this portion of the motion cues as shown in fig- 
ure 2. Appendix A contains plots of the motion base response 
in six axes of motion. 

The effective time delay associated with motion cueing 
logic was 89+16 milliseconds. This delay was determined 
in motion base tests (reference 3) and represents the average 
time from a discrete input to the host processor, command- 
ing a step input in the heave axis, until motion was detected 
by on-board accelerometers. It included a one and one-half 
simulation cycle frame time (at 32 ms per frame) and 41 ms 
of motion system response. The time delay did not include 
dynamics of the simulated aircraft. 

Cockpit layout 

The interior of the VMS cockpit is a generic four-en- 
gine, two-seat transport class aircraft as shown in figure 3. 
Six head-down color CRT displays are provided, along with 
analog indicators for altitude and airspeed. 

Two inceptors (control sticks) were provided: the left- 
side pilot station included a left-hand McFadden side-stick 
inceptor, and the right side pilot station included a center- 
stick inceptor. The two pilots shared a center aisle stand 
containing four throttles. 

The length of the side-stick from pivot point to middle 
of the grip was 6.5 inches. The length of the center-stick 
from pivot point to middle of the grip was 30.5 inches. 

Characteristics of these inceptors and throttles are given 
in tables 2 and 3 below. Force gradients are with respect to 
the middle of the stick grip. Bandwidth is estimated from 
the frequency that yielded 90 degree phase lag between com- 
manded and actual displacement in a calibration test. 
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Table 2. VMS Inceptor Static Characteristics 


Inceptor & 
Axis 

Test 

Amp. 

(in.) 

Band- 

width 

(rad/s) 

Break- 

out 

(lbs.) 

Force 

Gradient 

(lbs./in.) 

Neutral 

Position 

(degrees) 

Maximum 

Deflection 

(in.) 

Side Stick 
Pitch 

+0.16 

32 

1.4 

11 

9 fwd 

±1.4 

Side Stick 
Roll 

+0.15 

26 

1.4 

4 left 

5 right 

20 right 

±1.4 

Center 
Stick Pitch 

±0.20 

27 

2 

5.4 fwd 
6.8 aft 

0 

±3.6 

Center 
Stick Roll 

±0.40 

20 

2 

4.4 

0 

±3.6 

Pedals 

Yaw 

- 

- 

14 

65 

0 

±2.8 


Table 3. VMS Throttle Static Characteristics 


Throttle 
(numbered 
from left to 
right) 

Avg. 

Forward 

Friction 

(lbs.) 

Avg. Aft 
Friction 
(lbs.) 

Travel 

(inches) 

1 

3.7 

4.5 

8.6 

2 

4.9 

5.3 

8.6 

3 

5.6 

7.3 

8.6 

4 

3.7 

4.3 

8.6 





Figure 4. LaRC VMS Left-seat horizontal field- 
of-view. Right seat field-of-view same but 
reflected about centerline. 


Out-the-window displays 

Four mirror-beam-splitter-type displays were provided 
for out-the-window presentations, arranged with one display 
in front of each pilot, and one display on each side of the 
cockpit for a combined field-of-view for each pilot of ap- 
proximately 24 degrees vertical by 1 03 degrees horizontal 
(including a vertical seam between the front and side dis- 
play of 21 degrees) as shown in figure 4. These displays 
were based on XKD model 8000 monitors and each pre- 
sented an collimated image of 1023 pixels by 944 scan lines 
in front of the design eyepoint. 

Due to reuse of convex mirrors designed for another 
cab, the design eyepoint was located approximately 12 inches 
in front of the normal eyepoint as shown in figure 4. This 
results in minor magnification of the forward image, and a 
reduction in the exit pupil diameter - in effect, a portion of 
the forward image is masked by the outline of the mirror- 
beam-splitter at the normal eyepoint. 

Image generator 

The image generator used to provide out-the-window 
visual imagery was an Evans & Sutherland ESIG 3000 sys- 
tem. Terrain and culture data included a 200 nm by 200 nm 
area surrounding the Denver, Colorado International Airport. 
The database included geospecific photo texture within 20 
nm of the airport. For the puiposes of this research project, 
the altitude of the entire database was adjusted so that the 
threshold of the primary runway used (runway 35L) was set 
at mean sea level. The scene presented to the pilot was re- 
freshed at 60 Hz. 

The effective time delay associated with visual cueing 
was determined by testing (reference 3) and analysis to be 
108±24 milliseconds. This represents the amount of time 
required from a change in input at the cockpit, through the 
simulation host, running at a 32 ms frame time, to the comple- 
tion of the first visual field of the interlaced out-the-window 
display that reflected this change. The time delay did not 
include dynamics of the simulated aircraft. 

XVS symbology 

A representation of notional External Vision System 
(XVS) symbology was provided to the pilot. This represen- 
tation changed depending upon flight task. The basic sym- 
bology set is contained in figure 5 along with photographs 
of the actual symbology with and without out-the-window 
imagery. The out-the-window imagery was usually present 
during testing, with the exception of simulated Instrument 
Meteorological Conditions (IMC) when the ground imagery 
was obscured by simulated low visibility and clouds. 
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XVS mode annunciator 


TS 


TS 


S 


Tail scrape indicator 



(a) Basic symbology set of XVS 



(b) XVS symbology as implemented in LaRC.4 
(cruise conditions) 



(c) XVS symbology as implemented in LaRC.4 
(landing flare) 



(c) XVS symbology as implemented in LaRC.4 
(landing approach) 



(d) Pilots’ forward image with XVS symnbology 
(landing approach) 


Figure 5. Forward field-of-view with overlaid XVS symbology. 
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Heads-down displays 

Six 8” XYtron color cathode-ray tubes were mounted 
in the instrument panel of the VMS cockpit. These were used 
to present additional information to the flight crew, includ- 
ing a Primary Flight Display (figure 6), a Control Surface 
and Engine Display (figure 7) and a Navigation Display (fig- 
ure 8). If required for the task, a Velocity-Altitude Display 
(figure 9) was available. The Pilot Not Flying was provided 
with a Trim Display (figure 10) to verily the proper simula- 
tion trim state and that cockpit controls were in the proper 
position prior to beginning each simulation run. 


The math model of Reference H included aerodynam- 
ics, propulsion, landing gear, mass and inertia, sensor, and 
control law models. This model also included a representa- 
tion of structural dynamics of the aircraft, including the cou- 
pling between the aerodynamics, the actuators, and the elas- 
ticity of the vehicle (called dynamic aeroservoelastic, or 
DASE, effects). The first 40 modes of structural dynamics 
(the lowest 20 modes in frequency in both the symmetric 
and asymmetric axes) were modeled. These effects were 
modeled using a technique described as the “P-transform” 
method (see reference 4). The number of modes of each 
axis included in the simulation could be set at run time by 


Model description 

The mathematical description of the supersonic trans- 
port being investigated and much of the FORTRAN soft- 
ware used was produced by The Boeing Company. It was 
designated as Boeing Reference H, and was the fourth re- 
lease of the simulation model (Cycle 4). Modifications re- 
ceived after the initial release of the Cycle 4 model brought 
the Boeing-supplied software to a revision level designated 
as C. Additional software such as support routines, equa- 
tions of motion, and some elements of the structural mode 
control were written at NASA Langley. 





Figure 6. Primary Flight display. 



Figure 7. Control Surface and Engine display. 



cossc 

3S9 

/'a 

>> 


‘ 


— 25 — 






* .v y 

J 

$ 

is 372 


Figure 8. Navigation display. 



Figure 9. Velocity- Altitude display. 



Figure 10. Simulation trim display. 
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the console operator, including having no modes active (re- 
sulting in a quasi-steady aeroelastic model with no dynamic 
structural effects). 

Modifications of revision C 

Additional modifications were made to the Boeing-sup- 
plied code beyond those contained in revision C. These modi- 
fications are described in Appendix B. 

Control laws 

Control laws for both maneuvering and structural mode 
control (SMC) were provided by industry researchers in the 
form of higher-level software designs. These control system 
designs were refined during the course of LaRC.4. The ini- 
tial set of control laws were labeled version 0.23. Additional 
versions of the control laws were developed during LaRC.4, 
up to version 0.44. The result of these modifications led to 
the control laws used in a later experiment. Use of a higher- 
level software description of the control laws allowed for 
rapid modification in the computer-aided control design en- 
vironment and rapid incorporation of the changed design in 
the simulation. However, a characteristic of the code gen- 
erator used in this process led to excessive propagation de- 
lays for the higher-frequency SMC dynamic components, 
requiring hand-coding of the SMC portion of the longitudi- 
nal control system. 

Simulation schedule 

The schedule of pilot and engineer participation in the 
LaRC.4 study is shown in Appendix C. 

Session and run Logs 

A summary of sessions conducted in the LaRC.4 simu- 
lation study is given in Appendix D. The complete run log 
can be found in Appendix E. 

Significant results 

Lateral inceptor parametric study refined 

Development of a set of design requirements for a next- 
generation supersonic transport aircraft is one of the goals 
of the guidance and flight controls element of the HSR pro- 
gram. This study provided an opportunity to continue to 
develop these requirements for the lateral-motion character- 
istics of the pilot’s control inceptor, in this case, a center- 
stick configuration, hi an attempt to determine the desired 
relationship between the force gradient (pounds of force per 
degree of travel) of the lateral-axis control inceptor and the 
initial roll-axis motion response of the cockpit (p-dot), this 


relationship was varied parametrically and the pilot was asked 
to perform two different maneuvers: a supersonic small- 
amplitude heading change, a subsonic small-amplitude head- 
ing change, and an instrument approach to landing with dis- 
crete lateral gusts. 

The LaRC.4 study was used by the industry researchers 
to develop and refine these maneuvers and the matrix of para- 
metric variations, as well as to develop a task performance 
scoring algorithm for the small-amplitude heading changes. 

Rapid control law incorporation 

Since part of the mission of LaRC.4 was to provide a 
maneuver control law (non-DASE) and structural mode con- 
trol (SMC) law development opportunity, the use of a com- 
mercial -off-the-shelf computer-aided control law design pro- 
gram was supported. This control law design tool included a 
code generator that would convert the control design (repre- 
sented on the designer’s screen as a block diagram) into 
ANSI-C source, which was then compiled and linked into 
the real-time simulation program. This capability allowed 
the control designers to develop and refine candidate con- 
trol laws in a computer environment that was useful and fa- 
miliar, and allowed rapid incorporation of these control laws 
in the real-time simulation (on the order of 15 minutes). 

P-transform method validated 

LaRC.4 successfully validated the use of P-transform 
modeling of structural dynamics in a real-time simulation. 
Some initial errors in the implementation model were dis- 
covered and corrected in the course of LaRC.4 which im- 
proved the fidelity of the DASE P-transform formulation. 

SMC implementation discrepancy uncovered 

The auto-code process worked well for the rigid-body 
control law development and refinement portion of LaRC.4, 
but some difficulties arose during the SMC development ef- 
fort. In general, SMC laws tended to become unstable when 
operated in the real-time simulation even though they per- 
formed well in the computer-aided design program. Stabil- 
ity could be regained in the simulation by reducing the over- 
all SMC forward-path gain by a factor of five. 

Eventually, the cause of the disparity was determined to 
be inherent in the structure of the machine-generated con- 
trol law software. The SMC control laws were developed in 
the continuous (analog) domain. The nature of the gener- 
ated code was a set of subroutines ( in object-oriented termi- 
nology, methods) that allowed the simulation executive to 
obtain a set of state derivatives associated with the internal 
states of the control law for a given values of input variables 
and state variables. The executive would then integrate the 
control law state variables, using the obtained state deriva- 


6 



tives, to obtain updated state variables. Another machine- 
generated method was then called to obtain a new set of out- 
put variables associated with the set of input variables and 
updated state variables. 

Although this parallel integration method introduced 
consecutive one-cycle propagation delays for dynamic ele- 
ments (such as filters and lags) placed in series, the propa- 
gation delay was of little consequence for rigid-body con- 
trol laws in which (1) only a few low-order dynamic ele- 
ments were placed in series or (2) the iteration rate of the 
control laws was two orders of magnitude (or more) higher 
than the highest frequency of the simulated system. How- 
ever, in the case of structural mode control, the use of higher- 
frequency dynamic elements of higher-order, placed in se- 
ries, made this propagation delay sufficiently long that sig- 
nificant phase delay occurred and instabilities resulted. 

The proof of this delay as the cause of the instability 
was demonstrated by reducing the time step of the simula- 
tion model to much slower than real-time; when the simula- 
tion was run with a smaller time step, the SMC control law 
worked at the nominal gain settings. 

A workable solution to the problem was the manual con- 
version of a candidate longitudinal SMC control law into 
software, in which the state derivatives were integrated for 
each dynamic element before being passed along to the next 
dynamic element in the control law. 

Demonstration of structural inode control 

Following the correction to the SMC implementation, a 
successful set of structural mode control laws, working in 
concert with the rigid body control laws, was demonstrated 
with a pilot, with motion cueing on, for series of landing 
approaches. During this session, the SMC control laws were 
turned off, and then on, with significant improvement in both 
handling qualities and quality of ride noted by the research 
pilot. This was the first successful demonstration of SMC 
laws in the HSR program with a pilot-in-the-loop, 

Web- based access to results 

At the conclusion of each run, time-history data of se- 
lected parameters was stored on a computer disk in a digital 
format. This disk was shared by two computers - the real- 
time host processor, which wrote the data on the disk, and 
another computer that was running a World-Wide Web file 
server application. This Web server included the capability 
of providing secure, encrypted transfer of the time-history 
data using SSL protocols layered on top of the conventional 
TCP/IP Internet transfer protocols. 

By accessing a Web page on the secure server, an au- 
thorized remote user could browse through an up-to-date copy 


of the run log (also stored electronically at the end of each 
session) and select individual runs of interest, using a vari- 
ety of selection criteria. The runs of interest could then be 
downloaded securely over the Web to the user’s remote com- 
puter for further analysis. 

The selection criteria included: run number, pilot name, 
inceptor type, test type or flight card, control law configura- 
tion, motion cueing on or off, and DASE effects on or off. 

The data format used was the ASC2 format of Dry den’s 
GetData utility program (see reference 5). 

Another utility program written by one of the authors, 
tranth, was available for downloading as well. This program 
allowed the conversion (on UNIX computers) of ASC2 for- 
mat to other data formats (including Matlab binary and tab- 
separated values formats) for ease of import into various 
analysis programs. 

Remote participation in simulation 

LaRC.4 marked the first use of teleconferencing in LaRC 
simulations. During the evolution of this study, industry re- 
searchers requested follow-on simulator sessions after they 
had returned to their home offices on the West Coast. In one 
instance, an SMC designer wanted to get a pilot’s opinion of 
a change made in the SMC laws; in another instance, han- 
dling qualities engineers wanted to have a pilot assess a 
change made in the small-amplitude heading change maneu- 
ver. In both cases, the researchers were connected via tele- 
phone to the cockpit intercom. Although the visual scene 
was not presented to the remote researchers, they were able 
to converse with the pilot and obtain a description of the 
behavior and characteristics of the simulated aircraft ver- 
bally, as well as ask questions and clarifications of the pilot. 

hi all cases, time-history data from the simulation runs 
were available immediately after the event. 

This capability was useful in that it extended the evalu- 
ation time available to the researchers following the end of 
their on-site visit to Langley. 
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Appendix A: Langley VMS Motion Cueing System Response Plots 

The dynamic characteristics of the Langley Visual Motion Simulator were measured at the conclusion of LaRC.4 by 
driving the motion base, with and without dynamic compensation. To excite the motion platform within the linear response 
range of the system, a sine wave position command with the following amplitudes (+/— from zero ) was used as the input in 
each axis. Table A-l below gives the amplitude of excitation in each axis at each test frequency. 


Table A- 1 . Motion Response Input Sine Wave Amplitudes 



Frequency, rad/sec (Hz) 

Axis 

0.1 0.15 0.25 0.4 0.7 1.0 1.5 2.5 3.14 4.0 

(0.016) (0.024) (0.040) (0.064) (0.111) (0.16) (0.24) (0.40) (0.5) (0.64) 

6.28 

(1.0) 

9.0 

(1.43) 

12.6 

(2.0) 

Roll, 

Pitch, 

Yaw 

(deg) 

1.2 

1.0 

0.8 

0.5 

Surge, 

Sway, 

Heave 

(inch) 

1.0 

0.25 


The output of the motion base was measured from leg position sensors and converted back to equivalent motion base 
position. Figure A- 1 below shows a schematic representation of the test setup. 



Figure A-l. Method of obtaining motion base responses. 


The figures on the next three pages show the Bode response characteristics for this experiment for each axis of the 
VMS motion platform. The time delay depicted on each of the plots is derived by dividing the phase angle, in degrees, by 
the respective frequency, scaled to degrees per second. This is a separate measurement of the motion system time delay, 
which is frequency dependent and includes motion system dynamic characteristics, from the motion system transport delay 
measured in reference 3. 
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Appendix B: Modifications to Simulation Model subsequent to revision C 

The following discussion lists the changes applied to the Boeing-supplied implementation of the Cycle 4 Ref. H simulation 
model. Some of the changes are enhancements which support the use of the vehicle model for certain tasks (such as the 
zeroing of the engine inlet flow angularity at subsonic Mach numbers, which supports cross-wind takeoff tasks). Other 
changes are corrections to the implementation as originally delivered or subsequently revised (such as the actuator param- 
eter changes), or corrections to local interface routines which are necessary to work with Cycle 4 (such as the total pressure 
calculation for the engine model). 

This list was prepared by Gus Taylor of the Unisys Corp. 

Control and hinge moment models 

The Revision C versions of contrl mixer. f, flaps. f and ftc blk.f are used instead of the initially delivered versions. (SMR 98-15) 

Capabilities are added to the 3rd order actuator to bypass initialization in reset (bypass by default), to select the number of 
subframes (10 by default), and to select between Euler and AB2 integration (AB2 by default). Affects actuator 3rd. f. 
(SMRs 98-10 and 98-17) 

In the 3rd order actuator, position limits (LIM1 and LIM2) are applied to ACT POS, and the rate limit (rl) is applied to 
act_rate and cntl_dot. Affects actuator 3rd. f( SMR 98-10) 

The number of actuators used by ride control vanes is corrected (four rather than three). Affects actblkf (DR #42) 

Structural mode control command values are added to the rudder segment #1 command (CMND ( 4 ) ), the trailing edge flap 
segment #2 command (cmnd ( 7 ) ) and the trailing edge flap #7 command (cmnd ( 1 3 ) ), if the corresponding flags are true, 
before limiting. The values applied to the two flaps are negatives of one another. Affects contrl mixer.f . (cycle 4 control 
system) 

The sign of the limiting value of latrc_c is reversed: DTE7 c - clim(2, 9) rather than cl im (2, 9) - dte7_c. Af- 
fects contrl jnixerf. (DR #43) 

THE FOLLOWING CHANGES ARE NOT KNOWN TO BE NECESSARY FOR THE MODEL TO FUNCTION CORRECTLY 

Module contrl advan.f is not used. Control systems brought forward from cycle 3 are used instead. Affects cntrl execf 
(initial baseline) 

Initializations from hmoment blk.f are performed outside of the hinge moment model. The only delivered hinge moment 
module used in the real-time simulation is hmoment. f. 

The hinge moment model is called from subroutine ACT EXEC, and is called on every other pass only. Affects act exec. f. 
(initial baseline, SMR 98-10 and DR #40) 

Aerodynamic hinge moments (hmaer (i) ) are updated in subroutine actuator3 as well as in subroutine HMOMENT. Af- 
fects code in act exec f, actuator 3rd. f and hmoment. f. (SMR 98-10) 

Selection among multiple sets of 3rd order actuator parameters (fast, mid and baseline) is enabled. Affects act blk.f, 
act exec. f and actuator Srd.f (see attached code extracts). (SMR 98-12) 

Gearing is adjusted for the scaling of control inputs in the real-time program. Affects code in contrl iof. (initial baseline) 

Calculations of symmetric and antisymmetric components of control surface motions for use by the LaRC DASE model are 
added. (The LaRC DASE model was brought forward from cycle 3, and is distinct from the p-transform DASE model that 
was part of the delivered cycle 4 code.) Affects act exec. f. (initial baseline) 

Engine/propulsion model 

Updated versions of pro blk tab. f pro in buff, pro inlet cont.f and pro rnodel ./, received 2/17/98, are used in place of 
those initially delivered. 

The model requires that the total pressure (PT AR) input to it be calculated as follows regardless of Mach: 
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PT_AR = PRESS* (1.0 + 0 . 2*MACH_AR**2 ) **3 . 5 

where press is the ambient static pressure. This calculation produces the free-stream total pressure. Previously, the real- 
time simulation had used this calculation for the total pressure behind a normal shock. 

The percent net thrust (pcfnI through PCFN4) calculations are protected against division by zero. Affects engine out. f. 
(initial baseline) 

Variable thrstc O is not calculated in the engine model. Affects engine out. f. (initial baseline) 

Variables inic_PRO_R and comp_pro_r are used in place of imoden. Affects pro in buff, (initial baseline) 

Variables trm_thrust (i ) are not set as in the delivered code, but instead are set equal to the corresponding thrcmd_fcs ( i ) 
value on every pass in reset. Affects engine inf . (DR #7 and SMR 98-02) 

The low-speed engine inlet flow angularity calculation is modified by setting and applying gain XMWASH. Affects code in 
projniet_cont.fi see attached code extracts). (SMR 98-06) 

THE FOLLOWING CHANGE IS NOT KNOWN TO BE NECESSARY FOR THE MODEL TO FUNCTION CORRECTLY 

Calculations of the combined upwash and sidewash angles for each engine ( inlang ( i ) ), and of the limited inlet upwash/ 
sidewash tolerance for each engine (unsttol (i) ) are added. Affects code in pro inlet cont.f (see attached code extracts), 
(initial baseline and DR #4) 

P-transform DASE model 

Code from the Revision B version of ase Jlex.f and the Revision C version of aseinterp.f is used instead of that from the 
versions initially delivered, (revision b and SMR 98-15) 

The calculation of az_data ( i ) is reversed in sign, so that the sign convention of the az_data values is +down, consistent 
with the sign convention of azbi_ir and with the simulation data base document. Affects code in ase outf. (initial 
baseline) 

Signs are reversed in assigning gust velocities to U vector elements: 

U_S YM (22) = -WG 
and 

U_ASYM (26) = -VG 

This is necessary because the sign conventions assumed by the DASE model are +downward and +left-to-right. Affects 
ase inf (DR #25 and SMR 98-24) 

Inconsistent hardcoded input filename lengths are eliminated, and the capability is added to read 15 -character columns as 
well as 16-character columns from input files. Affects ase read ab f and ase read cf. (SMR 98-15) 

The gust velocity indices are corrected and parameterized. In the symmetric case, the index is NC_SYM* 3 + 1 = 22, and 
in the antisymmetric case it is NC_ASYM* 3 + 1 = 25. Affects code in ase Jlex.f. ( DR #30) 

In the DASE stabilizer doublet check case, it appears necessary to set longelv false, rather than true as hi the delivered 
dase Jong.dpt file, in order to match the delivered delevI plot (which is constant zero). 

THE FOLLOWING CHANGES ARE NOT KNOWN TO BE NECESSARY FOR THE MODEL TO FUNCTION CORRECTLY 

The capability is added to force point_ (a) sym to a selected value, regardless of the flight condition. When ride control 
vane and chin fin input data are in use, this permits a known set of input data to be selected regardless of how the data have 
been copied to fill in flight conditions outside of the ride control vane/chin fin domain. By default, the capability is turned 
on and the P0INT_ (A) sym values are forced to seven. Affects code in ase interp.f. (SMR 98-19) 

Actuator-provided values are used for control surface rates and accelerations, rather than computing rates and accelerations 
from deflections by backward difference. Gains defaulting to one are applied to the accelerations. Affects code in ase inf 
and ase Jlex.f. (SMRs 98-06 and 98-18) 

Gust acceleration terms are included in the DASE modal acceleration (zetad_ ( A) sym ( i ) ) calculations; gams defaulting 
to zero are applied to these accelerations. Affects code in ase inf and ase Jlex.f. (SMR 98-17) 
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The capability is added to inactivate (specifically, to residualize) all symmetric DASE modes above a selected number, and 
similarly for antisymmetric modes. (The capability to residualize modes by means of the nsel_ (a) sym ( i ) flags, con- 
tained in the delivered code, is retained as well.) Affects code in ase inf, rise Jlex.f and ase out f. (SMRs 98-09, 98-24 
and DR #28) 

Increments to mean-axis accelerations (ubid ir, vbid ir, etc.) due to DASE effects are calculated in trim in the same 
way as in operate, i.e., without using the abflex_ (a) sym matrices. In trim, the increments are not added into the accelera- 
tions, but the sum of their squares is added into the trim error. Similarly, DASE modal accelerations (ZETAD_ (A) sym (i ) ) 
are calculated in trim and the sum of their squares is added into the trim error. DASE modal deflections in trim are 
calculated directly from “A” and “B” matrices using equation 10.43 in the simulation data base document (p. 10-15), i.e., 
the abflex_ (A) sym matrices are no longer used for this purpose either. Affects code in ase exec. f, ase Jlex.f and ase interp. 
(SMR 98-23) 

Calculations of mean-axis load factors at the c.g. and of mean-axis external accelerations (NXBI_IR, NYBI_IR, nzbi ir, 
AXBI IR, AYBI_IR and AZBI IR) are performed outside of the DASE model. Affects code in ase out f (initial baseline) 

Gear/ground contact model 

Updating of stroke deflections (gp (i)) is inhibited if the zz element of the cosine matrix (zz_ER, limited as rzz in 
subroutine gear) is not positive. This avoids the calculation of large erroneous deflections when bank exceeds 90 deg. 
Affects gear./. (DR #34) 

THE FOLLOWING CHANGES ARE NOT KNOWN TO BE NECESSARY FOR THE MODEL TO FUNCTION CORRECTLY 

Many features are brought forward from earlier real-time cycles and incorporated into the gear model (see attached gear 
features). Affects code in gear.f, gear Jo. f and grd cntct.f. (initial baseline) 

Mass model 

The calculation of dycg_md is reversed in sign, so that the calculations using dycg_md in formomf, gear.f, and grdcntctf 
will be correct. Affects code in mass iof. (initial baseline) 

It is not known whether the fuel model is functional, because the model is not executed in the real-time simulation. No 
delivered code is affected; execution of the model is inhibited by setting variable PT_ANALYSIS to true. 

Aero model 

Variables CLTR and CMTR are looked up only once, from low speed/high speed tables. In the delivered code, they are 
looked up twice: in subroutine aero_tab (from low speed/high speed tables) and again in subroutine AERO_SPEC (from a 
special table). Affects aero spec. f. (initial baseline) 

Variables CLTGER_SZ and CLTGER_EZ are looked up from the CLTGER low-speed table (cltger ls.fda). In the delivered 
code, they are looked up from a special table (cltger. f da) which is identical to cltger ls.fda. The same is tiue of cmtger_S Z 
and CMTGER_EZ. Affects aero spec. f. (initial baseline) 

Equations-of-motion model 

The only executable module from this model used in the real-time simulation is formomf Code brought forward from 
earlier cycles is used in place of all other executable modules. Analysis indicates that the code in use is equivalent to the 
delivered model for the flat-earth case, i.e., the case for which the delivered model is valid. 

Atmosphere/gust model 

The delivered atmosphere/gust model is not used. 
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Code Extracts 

Listing D-l. Calculation and use of variable xmwash in pro_inlet_cont.fi 


C PROCESS pro_inlet_cont 1.4.1 

C * 

IF (sim_mode_pro_r . GE . 4) THEN 
C 

xmwash = (mach_ar - 1.0) / 0.2 
xmwash = min ( max (xmwash, 0 . 0) , 1.0) 

C * 


IF (inboardf_pro_l) THEN 
C * 

C PROCESS ic_local_ib 1.4. 1.3. 2 

C Inlet Local Conditions - Inboard 

C * 


angu_icm_pro_l (i_eng_pro_l) = (ql38 (i_eng_pro_l ) + 

& (ql3 9 ( i_eng_pro_l ) * angbaf_icf_pro_l ( i_eng_pro_l ) ) ) 

& * xmwash 


angs_icm_pro_l ( i_eng_pro_l ) = (ql35 (i_eng_pro_l) + 

& (ql37 (i_eng_pro_l) * angbaf_icf_pro_l (i_eng_pro_l) ) ) 

& * xmwash 

C * 

C END ic_local_ib 1 . 4 . 1 . 3 . 2 

C * 

END IF 

IF (outboardf_pro_l ) THEN 
C * 

C PROCESS ic_local_ob 1 . 4 . 1 . 3 . 3 

C Inlet Local Conditions - Outboard 

C * 


angu_icm_pro_l (i_eng_pro_l) = (ql25 (i_eng_pro_l) + 

& (ql26 ( i_eng_pro_l ) * angbaf_icf_pro_l (i_eng_pro_l) ) ) 

& * xmwash 


C 

C 

c 


angs_icm_pro_l (i_eng_pro_l) = (ql22 (i_eng_pro_l) + 

& (angbaf_icf_pro_l (i_eng_pro_l) * ql24 (i_eng_pro_l) ) ) 

& *xmwash 

it 

END ic_local_ob 1 . 4 . 1 . 3 . 3 

~k 

ENDIF 
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Listing D-2. Calculation of variables INLANG(i) and UNSTTOL(i) in pro_inlet_cont.fi 


k 

PROCESS ust_margin_2d 1.4. 1.9. 8 

Inlet Unstart Margin Calculation 

k 


unsttol (i_eng_pro_l) = ql21 (i_eng_pro_l) 

usmargex_imf_pro_l (i_eng_pro_l) = (qll7 (i_eng_pro_l ) . GT . 

& qll8 (i_eng_pro_l) ) .AND. ql 1 9 ( i_eng_pro_l ) 

* 

END ust_margin_2d 1 . 4 . 1 . 9 . 8 

-k 


k 

PROCESS in_ust_margin 1 . 4 . 1 . 4 . 8 

Inlet Unstart Margin Calculation 

-k 


unsttol (i_eng_pro_l) = q07 9 ( i_eng_pro_l) 

usmargex_imf_pro_l ( i_eng_pro_l) = (q075 (i_eng_pro_l ) . GT . 

& q07 6 ( i_eng_pro_l ) ) .AND. q077 (i_eng_pro_l ) 

k 

END in_ust_margin 1 . 4 . 1 . 4 . 8 

k 


inlang (i_eng_pro_l ) = qll7 (i_eng_pro_l) 

1000 CONTINUE ! end for loop 

k 

END pro_inlet_cont 1.4.1 

C * 


Listing D-3. Changes in act exec. f for multiple parameter sets: 


j = act_sel(i) 

call actuator3 (cntl (i) , cmnd (i) , cntl_dot (i) , cntl_ddot (i) , 

kk3 ( j , i ) , hmaer (i ) , hraact ( i) , chmd ( i) , chmdd(i), chmact(i), 
cchmact(i), hmpd (i ) , laap3 ( j , i) , 

nact 3 ( j , i ) , kvact3 ( j , i) , dim (l,i) ,clim(2,i) , rl3 ( j , i ) , psupply , 
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pact (i) , xv (i) , qv (i) , qf (i) , fl (i) , flp (i) , kp3 ( j , i) , 

c_t3 ( j , i ) , k_eqv3 ( j , i ) , k_str_i ( j , i) , i_hl3 ( j , i ) , blowdown , 

torque (i) , act_rate (i) , act_pos (i) ) 


Listing D-4. Changes in actuator 3rd. f for multiple parameter sets: 

SUBROUTINE ACTUATOR3 (CNTL, CNTL_CMD, CNTL_DOT , CNTL_DDOT , KK, HMAER, 

HMACT, CHMD, CHMDD, CHMACT , CCHMACT , HMPD , LAAP , 
NA, KV, LIM1, LIM2, RL, PS, PA, XVL, QV, QF, FL, FLP, 

KP , C T , K_EQV , K S TR_I , I_HL , BDOWN , TORQUE , 

ACT RATE, ACT POS) 


Thereafter: 

Input parameter laap is used in place of the product la*ap. 

The calculation of variable k_str is omitted, and input parameter k_str_i is used in its place. 


Listing D-5. Changes in actblkf for multiple parameter sets: 


DATA act sel / 24*1 / 


DATA rl3 / 25. 

25. 

, 

0. 

, 25. 

, 

! Stabilizer 


50. 

50 . 

, 

90. 

, 50. 

, 

! Left Hand Side Elevator 

#1 

50. 

50. 

/ 

90. 

, 50 . 

, 

! Right Hand Side Elevator #2 

50. 

50. 

/ 

90. 

, 50. 

/ 

! Rudder Panel #1 (lower) 


50. 

50. 

/ 

0. 

, 50. 

/ 

! Rudder Panel #2 (middle) 

50. 

50. 

/ 

0. 

, 50. 

/ 

! Rudder Panel #3 (upper) 


90. 

90 . 

/ 

90. 

, 90 . 

/ 

! Trailing edge flaperon 

#1 

90. 

90 . 

/ 

90. 

, 90. 

/ 

! Trailing edge flaperon 

#2 

90. 

90 . 

/ 

90. 

, 90. 

/ 

! Trailing edge flaperon 

#3 

15. 

15 . 


90. 

, 15. 

/ 

! Trailing edge flap #4 


15. 

15 . 


90 . 

, 15. 


! Trailing edge flap #5 


90. 

90 . 


90 . 

, 90 . 


! Trailing edge flaperon 

#6 

90. 

90 . 


90 . 

, 90 . 


! Trailing edge flaperon 

#7 

90. 

90 . 


90 . 

, 90 . 


! Trailing edge flaperon 

#8 

15. 

15 . 


0 . 

, 15 . 


! Leading edge flap #1 


15. 

15 . 


0 . 

, 15 . 


! Leading edge flap #2 


15. 

15 . 


0 . 

, 15 . 


! Leading edge flap #3 


15. 

15 . 


0 . 

, 15 . 


! Leading edge flap #4 


60. 

60 . 


0 . 

, 60 . 


! Spoiler Slot Deflector 

#1 

60. 

60 . 


0 . 

, 60 . 


! Spoiler Slot Deflector 

#2 

60. 

60 . 


0 . 

, 60 . 


! Spoiler Slot Deflector 

#3 

60. 

60 . 


0 . 

, 60 . 


! Spoiler Slot Deflector 

#4 

100 . 

90. 


90. 

, 90. 


! Canard 


100 . 

90. 


90. 

, 90. 

/ 

! Chin Vane 



DATA nact3/ 4 , 4 , 0 , 4 , 
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3,3,4, 

3,3,4, 

3.3.3, 

3.3.0, 

2 . 2 . 0 , 

2 , 2 , 2 , 

2 , 2 , 2 , 

2 , 2 , 2 , 

2.2.3, 

2.2.3, 

2 , 2 , 2 , 

2 , 2 , 2 , 

2 , 2 , 2 , 

6 , 6 , 0 , 

6 , 6 , 0 , 

6 , 6 , 0 , 

6 , 6 , 0 , 

2 , 2 , 0 , 

2 , 2 , 0 , 

2 , 2 , 0 , 

2 , 2 , 0 , 

4.4.3, 

2 , 2 , 2 , 


3 , 
3 , 
3 , 

3 , 
2 , 
2 , 
2 , 
2 , 
2 , 
2 , 
2 , 
2 , 
2 , 
6 , 
6 , 
6 , 
6 , 
2 , 
2 , 
2 , 
2 , 

4 , 
2 / 


DATA 


laap3/8 . 4220e+02 , 

1684.2 , 

0.0 , 

1684.2 , 

7 . 8952e+01, 

631.2 , 

103 .56, 

631.2 , 

7 . 8952e+01, 

631.2 , 

103.56, 

631.2 , 

8 . 6328e+01 , 

86.3 , 

24.34, 

86.3 , 

3 . 4 72 6e + 01 , 

34.7 , 

0.0 , 

34.7 , 

1 . 9971e+01, 

20.0 , 

0.0 , 

20.0 , 

1 . 6800e+01 , 

252.0 , 

85.35, 

252.0 , 

1 . 4001e+02, 

560.0 , 

172.9 , 

560.0 , 

2 ,0009e+02, 1200.0 , 

292.5 , 

1200.0 , 

7 . 4520e+01, 

74.5 , 

342.7 , 

74.5, 

7 . 4520e+01, 

74.5, 

342.7 , 

74.5, 

2. 0009e+02, 1200.0 , 

292.5 , 

1200.0 , 

1 . 4001e+02, 

560.0 , 

172.9 , 

560.0 , 

1 . 6800e+01 , 

252.0 , 

85.35, 

252.0 , 

5 . 7933e+01, 

57 . 93, 

0.0 , 

57.93, 

3 . 8640e+01 , 

38.64, 

0.0 , 

38.64, 

3 . 8 64 0e+01 , 

38 . 64, 

0.0 , 

38.64, 

5 . 7 933e+01 , 

57 . 93, 

0.0 , 

57.93, 

1 . 5200e+01 , 

15.2 , 

0.0 , 

15.2 , 

7 . 2000e+00, 

7.2 , 

0.0 , 

7.2 , 

7 . 2000e+00, 

7.2 , 

0.0 , 

7.2 , 

1 . 5200e+01 , 

15.2 , 

0.0 , 

15.2 , 

1 . 0000e+00, 

1301 . 3, 

1301 . 3, 

1301.3 , 

1 . 0000e+00, 

90.3, 

90.3, 

90.3 / 


DATA k str i 


/ 2 . 0956e-09 , 

2 . 24e-19, 

o 

o 

2 . 24e-19 

5. 0828e-08, 

51 . 62e-19, 

10. 91e-l 9 , 

51 . 62e-l 9 

5. 0828e-08, 

51 . 62e-19, 

10 . 91e-l 9 , 

51 . 62e-l 9 

2 . 8730e-07 , 

287 . 3e-19, 

46. 58e-19, 

287 . 3e-l 9 

5. 0326e-07, 509 ,29e-19. 

o 

o 

509.2 9e-l 9 

1 . 01 99e-06 , 1019. 9e-19, 

o 

o 

1019 . 9e-19 

1 . 27 69e-07 , 127 . 32e-19, 

13 . 29e-19, 

127 . 32e-l 9 

6. 4354e-08, 

64 . 35e-l 9 , 

6 . 7e-l 9 , 

64 . 35e-l 9 

3. 6544e-08 , 

36 . 38e-l 9 , 

3 . 82e-l 9 , 

36. 38e-19 

2 . 1417e-08, 

2 . 1417e-8, 

3. 36e-19, 

o 

o 

2 . 1417e-08, 

2 . 1417e-8, 

3. 36e-19, 

o 

o 

3. 6544e-08, 

36 . 38e-19, 

3 . 82e-19, 

36 . 38e-19 

6 . 4354e-08, 

64 . 35e-l 9 , 

6 . 7e-19, 

64 . 35e-19 
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u u 


c 


loop gain 
DATA 


1 . 

27 69e-07 , 

127 . 32e- 

19, 13 . 29e-19, 

127 

1 . 

4 0 64e-07 , 

1 . 408e 

-7, 

0.0, 

1 

2 . 

41 30e-07 , 

2 . 413e 

-7, 

0.0, 

2 

2 . 

41 30e-07 , 

2 . 413e 

-7, 

0.0, 

2 

1 . 

4064e-07 , 

1 . 408e 

-7, 

0.0, 

1 

7 . 

71 66e-06 , 

7 . 717e 

-6, 

0.0, 

7 

7 . 

71 66e-06 , 

7 . 717e 

-6, 

0.0, 

7 

7 . 

71 66e-06 , 

7 . 717e 

-6, 

0.0, 

7 

7 . 

71 66e-06 , 

7 . 717e 

-6, 

0.0, 

7 

7 . 

6394e-01 , 

8 . 73e- 

19, 

0.0, 

8 

7 . 

6394e-01 , 

12 . 73e- 

19, 

0.0, 

12 

s for 

3'rd order act. 

(in /deg) 



0176 

, .0176 , 

.0000 , 

.0176 , 



0088 

, .0088 , 

. 00489, 

.0088 , 



0088 

, .0088 , 

. 00489, 

.0088 , 



0088 

, .0088 , 

. 00489, 

.0088 , 



0088 

, .0088 , 

.0000 , 

.0088 , 



0088 

, .0088 , 

.0000 , 

.0088 , 



00489 

, .00489, 

. 00489, 

. 00489, 



00489 

, .00489, 

. 00489, 

.00489, 



00489 

, .00489, 

.00489, 

.00489, 



02936 

, .02936, 

.00489, 

.02936, 



02936 

, .02936, 

.00489, 

.02936, 



00489 

, .00489, 

.00489, 

. 00489, 



00489 

, .00489, 

. 00489, 

.00489, 



00489 

, .00489, 

.00489, 

. 00489, 



02936 

, .02936, 

. 00000, 

.02936, 



02936 

, .02936, 

. 00000, 

.02936, 



02936 

, .02936, 

.00000, 

. 02936, 



02936 

, .02936, 

.00000, 

.02936, 



00735 

, .00735, 

.00000, 

. 00735, 



00735 

, .00735, 

. 00000, 

. 00735, 



00735 

, .00735, 

. 00000, 

.00735, 



00735 

, .00735, 

. 00000, 

.00735, 



00000 

, .00489, 

. 00489, 

.00489, 



00000 

, .00489, 

.00489, 

.00489/ 




DATA kp3 / .00078, 

.00010 

, 

. 00000 


.00010 , 

. 00156, 

. 00010 

, 

. 00001 


.00010 , 

. 00156, 

.00010 

, 

. 00001 


.00010 , 

.0097 , 

.00010 

, 

. 00001 


.00010 , 

.0068 , 

.00010 

/ 

. 0000 


.00010 , 

.0065 , 

. 00010 

/ 

. 0000 


.00010 , 

. 00068, 

. 00010 

/ 

. 0000 


.00010 , 

.0028 , 

. 00010 

/ 

. 0000 


.00010 , 

.0023 , 

. 00010 


. 0000 


.00010 , 

.0000 , 

. 0000 


. 0000 


.0000 , 

.0000 , 

. 0000 


. 0000 


.0000 , 

.0023 , 

. 00010 


. 0000 


.00010 , 

.0028 , 

. 00010 


. 0000 


.00010 , 

. 00068, 

. 00010 


. 0000 


.00010 , 

.0000 , 

. 0000 


. 0000 


.0000 , 

.0000 , 

. 0000 


. 0000 


.0000 , 

.0000 , 

. 0000 


. 0000 


.0000 , 

.0000 , 

. 0000 


. 0000 


.0000 , 

.0000 , 

. 0000 


. 0000 


.0000 , 

.0000 , 

. 0000 


. 0000 


.0000 , 

.0000 , 

. 0000 


. 0000 


.0000 , 

.0000 , 

. 0000 


. 0000 


.0000 , 

.0000 , 

. 00001 


. 00001 


.00001 , 

.0000 , 

. 00001 


. 00001 


.00001 / 


32e-19 

408e-7 

413e-7 

413e-7 

408e-7 

717e-6 

717e-6 

717e-6 

717e-6 

73e-19 

73e-19 
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DATA 


DATA 


DATA 


c t3 / 


0 . 125E+05 
0 . 450E+03 
0 . 450E+03 
0 . 800E+02 
0 . 460E+02 
0 . 190E+02 
0 . 155E+03 
0 . 290E+03 
0 . 500E+03 
0 . OOOE+OO 
0 . OOOE+OO 
0 . 500E+03 
0 . 2 90E+03 
0 . 155E+03 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 .OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 .OOOE+OO 


0 . 209E+06 
0 . 803E+04 
0 . 803E+04 
0 . 142E+04 
0 . 803E+03 
0 . 332E+03 
0 . 271E+04 
0 . 524E+04 
0 . 916E+04 
0 . OOOE+OO 
0 . OOOE+OO 
0 . 916E+04 
0 . 524E+04 
0 . 271E+04 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . 1 92E+05 
0 . 873E+04 


0 . 218E+03 
0 . 204E+05 
0 . 204E+05 
0 . 361E+04 
0 . OOOE+OO 
0 . OOOE+OO 
0 . 838E+04 
0 . 166E+05 
0 . 290E+05 
0 . 497E+05 
0 . 497E+05 
0 . 290E+05 
0 . 166E+05 
0 . 838E+04 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . 1 92E+05 
0 . 873E+04 


0 . 209E+06 
0 . 803E+04 
0 . 803E+04 
0 . 142E+04 
0 . 803E+03 
0 . 332E+03 
0 . 271E+04 
0 . 524E+04 
0 . 916E+04 
0 . OOOE+OO 
0 . OOOE+OO 
0 . 916E+04 
0 . 524E+04 
0 . 271E+04 
0 . OOOE+OO 
0 .OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . OOOE+OO 
0 . 192E+05 
0 . 873E+04 


kvact3/ 


236.16 , 
44.26 , 
44.26 , 

48.41 , 
19.46 , 

11.22 , 
16.96 , 
141.32 , 
201.89 , 
12.54 , 
12.54 , 
201.89 , 
141.32 , 
16.96 , 

9.76 , 

6.51 , 

6.51 , 

9.76 , 

10.23 , 

4.85 , 

4.85 , 

10.23 , 
.0000 , 
.0000 , 


500.0 , 

360.0 , 

360.0 , 

48.41 , 

20.00 , 
11.22 , 

262.5 , 

580.0 , 

1230.0 , 

12.53 , 

12.53 , 

1230.0 , 

580.0 , 

262.5 , 

9.76 , 

6.51 , 

6.51 , 

9.76 , 

10.23 , 

4.85 , 

4.85 , 

10.23 , 

390.0 , 
268.75, 


0 . , 

307.5 , 

307.5 , 
72.60 , 
0 . , 
0 . , 
247.65 , 
500.50 , 
843.75 , 

989.00 , 

989.00 , 
843.75 , 
500.50 , 
247.65 , 
0 . , 
0 . , 
0 . , 
0 . , 
0 . , 
0 . , 
0 . , 
0 . , 

390.0 , 
268 .75, 


500.0 , 

360.0 , 

360.0 , 

48.41 , 

20.00 , 
11.22 , 

262.5 , 

580.0 , 

1230.0 , 

12.53 , 

12.53 , 

1230.0 , 

580.0 , 

262.5 , 

9.76 , 

6.51 , 

6.51 , 

9.76 , 

10.23 , 

4.85 , 

4.85 , 

10.23 , 

390.0 , 
268.75/ 


k_eqv3/ 


36. 4550e7 

34 . Ie7 

1 . 5030e7 

1 . 48e7 

1 . 5030e7 

1 . 48e7 

. 2659e7 

. 2659e7 

. 1518e7 

. 15e7 

. 0749e7 

. 0749e7 

. 5983e7 

. 60e7 

1 . 1871e7 

1 . 1871e7 

2 . 0905e7 

2 . 1000e7 

3 . 5670e7 

3 . 5670e7 


0 . , 

34 . Ie7 

7 . 0e7 , 

1 . 48e7 

7 . 0e7 , 

1 . 48e7 

1 . 64e7 , 

. 2 65 9e7 

0 . , 

. 15e7 

0 . , 

. 0749e7 

5 . 7 5e7 , 

. 6 0e7 

11 . 4e7 , 

1 . 1871e7 

2 0 . Oe7 , 

2 . 1000e7 

22 . 75e7 , 

3 . 5670e7 



3 . 5670e7 , 
2 . 0905e7 , 
1 . 1871e7 , 
. 5983e7 , 
. 5432e7 , 
. 3166e7 , 
. 3166e7 , 
. 5432e7 , 
. 0099e7 , 
. 0099e7 , 
. 009 9e7 , 
. 0099e7 , 
.0000 
.0000 


3 . 5 67 0e7 , 
2 . 1 000e7 , 
1 . 1871e7 , 
. 60e7 , 
. 5432e7 , 
. 31 66e7 , 
. 31 66e7 , 
. 5432e7 , 
. 0099e7 , 
. 0099e7 , 
. 0099e7 , 
. 0099e7 , 
8 . 75e7 , 
6 . 0e7 , 


22.7 5e7 , 
20 . 0e7 , 
11 . 4e7 , 
5 . 7 5e7 , 
0 . , 
0 . , 
0 . , 
0 . , 
0 . , 
0 . , 
0 . , 
0 . , 
8 . 75e7 , 
6 . 0e7 , 


3 . 5670e7 , 
2 . 1000e7 , 
1 . 1871e7 , 

. 60e7 , 

. 5432e7 , 

. 31 66e7 , 

. 31 66e7 , 

. 5432e7 , 

. 0099e7 , 

. 0099e7 , 

. 0099e7 , 

. 0099e7 , 
8 . 75e7 , 
6 . 0e7 / 


DATA i hl3/ 35644. 24e3 

, 35644. 20e3 , 

0 . , 

35644. 24e3 , 

1469. 52e3 

, 1469. 52e3 , 

1469. 52e3 , 

1469. 52e3 , 

1469. 52e3 

, 1469. 52e3 , 

1469. 52e3 , 

1469. 52e3 , 

259. 98e3 

, 259 . 98e3 , 

259 . 98e3 , 

259 . 98e3 , 

148 . 38e3 

, 148 . 38e3 , 

0 . , 

148 . 38e3 , 

7 2 . 32e3 

, 72 . 32e3 , 

0 . , 

72 . 32e3 , 

585. 02e3 

, 585 . 02e3 , 

585 . 02e3 , 

585 . 02e3 , 

1160. 65e3 

, 1160. 65e3 , 

1160. 65e3 , 

1160. 65e3 , 

2044 . 02e3 

, 2044. 02e3 , 

2044. 02e3 , 

2044. 02e3 , 

3487 . 60e3 

, 3487. 60e3 , 

3487. 60e3 , 

3487. 60e3 , 

3487 . 60e3 

, 3487. 60e3 , 

3487. 60e3 , 

3487. 60e3 , 

2044 . 02e3 

, 2044. 02e3 , 

2044. 02e3 , 

2044. 02e3 , 

1160. 65e3 

, 1160. 65e3 , 

1160. 65e3 , 

1160. 65e3 , 

585 . 02e3 

, 585 . 02e3 , 

585 . 02e3 , 

585 . 02e3 , 

531 .15e3 

, 531.1 5e3 , 

0 . , 

531 . 15e3 , 

309 . 53e3 

, 309 . 53e3 , 

0 . , 

309 . 53e3 , 

309 . 53e3 

, 309 . 53e3 , 

0 . , 

309 . 53e3 , 

531 . 1 5e3 

, 531 . 15e3 , 

0 . , 

531 . 15e3 , 

9. 67e3 

, 9 . 67e3 , 

0 . , 

9 . 67e3 , 

9. 67e3 

, 9 . 67e3 , 

0 . , 

9 . 67e3 , 

9. 67e3 

, 9 . 67e3 , 

0 . , 

9 . 67e3 , 

9 . 67e3 

, 9 . 67e3 , 

0 . , 

9 . 67e3 , 

. 0000 

, 1380. 00e3 , 

1380. 00e3 , 

1380. 00e3 , 

. 0000 

, 630 . 00e3 , 

630 . 00e3 , 

630 . 00e3 / 


Cycle 3 features brought forward into the Cycle 4 gear/ground contact model 


Antiskid braking 
Runway roughness model 

Capability to select between different strut spring models 

Calculation of gear altitudes 

Calculation of projected tailscrape pitch angle 

Nose wheel steering features: 

A tiller capability is provided. 

The gain strp is found by interpolation in reset just as in operate, i.e.. there is no special reset value 

Variable rudpedp cp is assumed to be already normalized, and no trim value is provided for, i.e., vari- 
ables ped gain and dg tr are not used. 


Brake features: 


Variables rbped and lbped are set either to the cab input or to the trim value, not to the sum. and are 


22 



assumed to be already normalized. 

rbped, lbped and bped are scaled rather then being limited. 

Side force computation features: 

deltau and deltay ate calculated (not zeroized) in reset. 

Past values: 

fygg ( i ) values are zeroized in reset. 

pug (i) , pvg (i) andFXGG(i) values are initialized in reset only if on_ground_tr is false, in which 
case they are set to the last values held when on ground tr was true. 


Initializations differing from Cycle 4: 


AUTBDZ 

= 

4 

AMUBIC 

= 

0.016 

AUTBRK 

= 

. FALSE 

BRAKE GEARING (1) 

= 

2 

FLAG 

= 

0 

GDEF (1) 

= 

0.0 

GNBPED 

= 

0.0 

NWSG2 

= 

15.0 

NWSGS2 

= 

1000.0 

NWSRLM 

= 

60.0 

RLGU 

= 

9.0 


autobrake deceleration rate = 1.8 (DECFS (DECEL) in Boeing code) 
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Appendix C: Simulation Schedule 


Wk 

Dates 

(1998) 

Purpose 

Engineers 

Pilots 

0 

3/30 - 4/3 

LaRC checkout 

Jackson 

Rivers 

1 

4/6 - 4/10 

FCS checkout 

Rowhani 

Williams 

Rivers 

Norman 

2 

4/13 - 4/17 

Kraft 

3 

4/10 - 4/24 

SMC checkout 

Borland 

Rowhani 

Goldthorpe 


■ 

4/27 - 5/1 

DASE verification. 
Roll Sensitivity 
preparations 

Coleman 

Hardy 

5 

5/4 - 5/8 

von Klein 


Rivers 

6 

5/11-5/15 

DASE validation 



7 

5/18 - 5/22 

Roll Sensitivity Study 

Rossitto 

Shweyk 


8 

5/26 - 5/29 

Borland 

Verstynen 

9 

6/1 - 6/5 

SMC checkout 

Goldthorpe 

Nagaraja 

Rowhani 

Barrett 

Norman 
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Appendix D: Session Log 


Session 

Number 

Run 

Date 

Session 

Start 

Session 

Finish 

Elapsed 

Time 

Start 

Run 

End 

Run 

Run 

Count 

Engineer 
or Pilot 

Notes 

1 

4/8 

8:30 

11:00 

2:30 

76 

111 

36 

Rowhani 


2 

4/8 

14:26 

14:54 

0:28 

113 

116 

4 

Williams 


3 

4/8 

16:29 

17:00 

0:31 

120 

130 

11 

Rowhani 


4 

4/9 

8:24 

10:17 

1:53 

136 

172 

37 

Rowhani 


5 

4/9 

12:11 

12:52 

0:41 

178 

191 

14 

Norman 


6 

4/9 

13:41 

14:00 

0:19 

193 

202 

10 

Williams 


7 

4/9 

16:15 

16:26 

0:11 

213 

218 

6 

Williams 


8 

4/10 

8:23 

9:01 

0:38 

224 

233 

10 

Rowhani 


9 

4/10 

9:01 

10:05 

1:04 

234 

253 

20 

Norman 


10 

4/10 

10:08 

10:10 

0:02 

254 

255 

2 

Raney 


11 

4/10 

10:11 

11:00 

0:49 

256 

260 

5 

Williams 


12 

4/10 

13:22 

14:26 

1:04 

263 

277 

15 

Rivers 


13 

4/10 

14:47 

16:07 

1:20 

278 

286 

9 

Rivers 


14 

4/13 

9:10 

10:00 

0:50 

301 

307 

7 

Kraft 


15 

4/13 

15:31 

16:02 

0:31 

344 

356 

13 

Williams 


16 

4/13 

16:16 

16:55 

0:39 

357 

368 

12 

Kraft 


17 

4/14 

9:17 

10:06 

0:49 

371 

388 

18 

Kraft 


18 

4/14 

10:41 

10:54 

0:13 

389 

393 

5 

Williams 


19 

4/14 

14:35 

15:24 

0:49 

396 

409 

14 

Kraft 


20 

4/14 

16:00 

16:22 

0:21 

410 

415 

6 

Williams 


21 

4/15 

10:43 

11:00 

0:17 

432 

436 

5 

Rivers 


22 

4/15 

14:30 

15:22 

0:52 

440 

450 

11 

Rivers 


23 

4/16 

10:27 

11:06 

0:39 

470 

480 

11 

Rivers 


24 

4/16 

15:07 

16:03 

0:56 

489 

501 

13 

Rivers 


25 

4/17 

9:28 

11:04 

1:36 

507 

537 

31 

Rivers 


26 

4/17 

14:34 

14:40 

0:06 

539 

540 

2 

Kraft 


27 

4/17 

14:48 

15:43 

0:55 

541 

556 

16 

Norman 


28 

4/21 

9:04 

9:35 

0:31 

604 

621 

18 

Borland 


29 

4/21 

15:18 

17:11 

1:53 

625 

635 

11 

Borland 


30 

4/22 

10:56 

11:30 

0:34 

673 

676 

4 

Goldthorpe 


31 

4/22 

12:59 

14:19 

1:20 

679 

697 

19 

Elenderson 


32 

4/22 

14:33 

14:40 

0:07 

698 

701 

4 

Borland 


33 

4/22 

15:28 

15:32 

0:04 

708 

710 

3 

Elenderson 

got lat/dir SMC working 

34 

4/22 

16:14 

16:43 

0:29 

718 

739 

22 

Borland 

got long SMC working 

35 

4/22 

16:45 

17:00 

0:14 

730 

746 

17 

Elenderson 

removed notch and 
bandpass filters; both axes 
under SMC 

36 

4/23 

8:40 

8:59 

0:18 

767 

774 

8 

Elenderson 


37 

4/23 

9:38 

10:12 

0:34 

781 

793 

13 

Rowhani 

potential gamma checkout 

38 

4/23 

10:14 

11:00 

0:45 

794 

804 

11 

Elenderson 

Lat/Dir SMC 

39 

4/23 

15:34 

15:46 

0:12 

816 

819 

4 

Goldthorpe 

Lat/Dir SMC, ECS vU.36 

40 

4/23 

15:46 

16:40 

0:53 

820 

839 

20 

Borland 

Checkcases 

41 

4/27 

10:14 

11:30 

1:16 

874 

875 

2 

Coleman 

prot. 

42 

4/27 

14:39 

16:20 

1:40 

879 

904 

26 

Coleman 

env. prot; rolls 

43 

4/29 

16:54 

16:59 

0:05 

938 

941 

4 

Elardy 

Flare cue tuning 

44 

4/30 

10:53 

11:05 

0:11 

949 

949 

1 

Elardy 

Gamma tuning; roll 
snubbing 

45 

4/30 

14:10 

16:04 

1:54 

954 

978 

25 

Borland 

SMC checkout motion 

46 

4/30 

16:16 

16:50 

0:34 

979 

986 

8 

Hardy 

DASE demo with Gordon 

47 

4/30 

16:51 

16:58 

0:06 

987 

988 

2 

Derry 

Demo to Julie Mikula 

48 

5/1 

9:00 

9:24 

0:24 

996 

1008 

13 

Borland 

New Bl_. 13 checkout 

49 

5/1 

9:30 

9:50 

0:20 

1009 

1013 

5 

Bailey 

Center stick sweeps 

50 

5/1 

9:58 

10:59 

1:00 

1014 

1039 

26 

Hardy 

Flare cue tuning for Ames. 

51 

5/1 

14:34 

15:02 

0:28 

1043 

1051 

9 

Borland 

SMC checkout 

52 

5/1 

15:34 

16:34 

1:00 

1054 

1081 

28 

Hardy 

Flare bias, and A/T 
transients 

53 

5/4 

9:00 

10:45 

1:45 

1087 

1097 

11 

von Klein 

roll hold checkout 

54 

5/5 

10:07 

10:59 

0:52 

1112 

1137 

26 

von Klein 

roll hold checkout 
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Session Log (cont’d) 


Session 

Run 

Session 

Session 

Elapsed 

Start 

End 

Run 

Engineer 


Number 

Date 

Start 

Finish 

Time 

Run 

Run 

Count 

or Pilot 

Notes 

55 

5/6 

8:45 

9:50 

1:05 

1140 

1155 

16 

von Klein 

roll control laws 

56 

5/4 

9:00 

10:45 

1:45 

1087 

1097 

11 

von Klein 

roll hold checkout 

57 

5/5 

10:07 

10:59 

0:52 

1112 

1137 

26 

von Klein 

roll hold checkout 

58 

5/6 

8:45 

9:50 

1:05 

1140 

1155 

16 

von Klein 

roll control laws 

59 

5/7 

8:30 

11:00 

2:30 

1168 

1186 

19 

von Klein 

roll control laws 

60 

5/8 

8:30 

10:49 

2:19 

1195 

1219 

25 

Norman 

checkout of roll gains and 
stick sensitivity 

61 

5/14 

14:45 

15:28 

0:43 

1267 

1275 

9 

Norman 

check rudder RL; used 
conference call with T. 










vonK 

62 

5/18 

10:00 

11:00 

1:00 

1293 

1296 

4 

Rossito 

sweep of controls 

63 

5/18 

15:35 

17:00 

1:24 

1296 

1322 

27 

Shweyk 


64 

5/19 

9:40 

10:42 

1:02 

1330 

1341 

12 

Rossitto 


65 

5/20 

9:28 

11:02 

1:34 

1356 

1379 

24 

Norman 


66 

5/20 

14:52 

15:19 

0:27 

1385 

1389 

5 

Rossitto 


67 

5/21 

10:30 

11:00 

0:30 

1402 

1405 

4 

Rossitto 


68 

5/21 

14:42 

17:00 

2:18 

1408 

1436 

29 

Norman 

begin roll cards 

69 

5/22 

9:16 

11:02 

1:46 

1439 

1468 

30 

Norman 

finished landing card 

70 

5/22 

14:46 

16:50 

2:03 

1477 

1505 

29 

Norman 

supersonic heading cards 

71 

5/26 

10:22 

11:00 

0:38 

1512 

1518 

7 

Rivers 

begin roll cards 

72 

5/26 

15:19 

17:03 

1:44 

1522 

1538 

17 

Rivers 

continued landing cards 

73 

5/27 

9:07 

11:01 

1:54 

1541 

1560 

20 

Verstynen 

subsonic roll 

74 

5/27 

15:17 

17:15 

1:57 

1564 

1583 

20 

Rivers 

finish landing cards 

75 

5/28 

8:55 

11:02 

2:07 

1587 

1600 

14 

Rivers 

supersonic roll cards 

76 

5/28 

15:00 

16:25 

1:24 

1603 

1617 

15 

Verstynen 


77 

5/28 

16:28 

17:05 

0:36 

1618 

1624 

7 

Rivers 

finish supersonic cards 

78 

5/29 

8:47 

8:55 

0:08 

1628 

1636 

9 

Otto 

roll check cases 

79 

6/1 

8:58 

11:30 

2:32 

1650 

1695 

46 

Goldthorpe 

rudder pulses 

80 

6/1 

12:52 

16:16 

3:24 

1700 

1730 

31 

Goldthorpe 

VMS checkout of lat/dir 
SMC 

81 

6/2 

8:30 

10:00 

1:30 

1731 

1748 

18 

Goldthorpe 

Gain margin investigations 

82 

6/2 

9:44 

11:00 

1:16 

1751 

1785 

35 

Borland 

long smc checkout 

83 

6/2 

11:15 

14:41 

3:26 

1788 

1830 

43 

Norman 

Flew with Steve G with 
lat/dir on motion 

84 

6/2 

14:52 

16:04 

1:12 

1831 

1847 

17 

Borland 

long smc checkout 

85 

6/3 

8:30 

8:54 

0:24 

1866 

1870 

5 

Norman 

examine SMC a2r 

86 

6/3 

10:57 

16:09 

5:12 

1870 

1901 

32 

Borland 

slowed down sim 4X 

87 

6/3 

16:27 

17:03 

0:35 

1901 

1906 

6 

Barrett 


88 

6/4 

8:30 

9:48 

1:18 

1921 

1935 

15 

Gregory 

first flight 

89 

6/4 

9:57 

10:40 

0:43 

1936 

1960 

25 

Borland 

George Hunt's SMC coding 

90 

6/4 

10:40 

11:07 

0:27 

1961 

1966 

6 

Barrett 

Tracked down duplicate 
flags 

91 

6/4 

14:31 

15:45 

1:13 

1969 

1990 

22 

Barrett 

showed numerical 
instability 

92 

6/4 

16:05 

17:30 

1:25 

1996 

2017 

22 

Borland 

invest, various param 
changes 

93 

6/5 

8:35 

9:25 

0:49 

2021 

2037 

17 

Barrett 


94 

6/5 

9:32 

11:02 

1:30 

2039 

2053 

15 

Borland 


95 

6/5 

14:44 

15:20 

0:35 

2056 

2059 

4 

Barrett 

flew with Mike N 

96 

6/5 

15:38 

16:25 

0:47 

2061 

2074 

14 

Gregory 


97 

6/8 

14:30 

15:44 

1:14 

2079 

2103 

25 

Norman 

Looked at hdg task changes 

98 

6/8 

16:06 

16:19 

0:13 

2109 

2115 

7 

Gregory 


99 

6/17 

14:30 

16:20 

1:50 

2144 

2167 

24 

Norman 


100 

6/17 

16:25 

17:00 

0:35 

2168 

2175 

8 

Rivers 


101 

6/24 

9:15 

9:30 

0:15 

2179 

2187 

9 

Otto 


Totals 

106.23 

hours 


1524 

runs 
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Appendix F: Task Identifiers and Flight Cards 

Each of the tasks listed in the run log (appendix E) corresponded to an initial conditions description file which con- 
tained initial values for simulator states, such as airspeed, altitude, heading, and other simulation model parameters. Table 
F-l below lists some of the initial conditions for these tasks, including a short identifier for the task. 


Table F-l. Task Names and Initial Conditions 


Task 

ID 

Task Name 

Weight 

(lbs) 

C.G. 

(%) 

Velocity 

Altitude 

(ft) 

2010 

Standard Acoustic Takeoff 

649,914 

48.1 

0 

0 

2011 

Alternate Acoustic Takeoff 

649,914 

48.1 

0 

0 

2030 

Programmed Lapse Rate Takeoff with Gusts 

649,914 

48.1 

0 

0 

3030 

Profile Climb 

649,914 

48.1 

0 

0 

3076 

Subsonic Deceleration 

384,862 

53.2 

345 kt 

15,000 

3084 

Heading Change in Late Supersonic Cruise 

384,862 

53.2 

M 2.4 

63,700 

3086 

Heading Change in Low Altitude Cruise 

384,862 

53.2 

345 kt 

15,000 

3384 

Heading Change in Supersonic Cruise - Small Inputs 

384,862 

53.2 

M 2.4 

64,700 

3385 

Heading Change in Subsonic Cruise - Small Inputs 

384,862 

53.2 

250 kt 

25,000 

4020 

Nominal Approach & Landing 

384,862 

53.2 

190 kt 

1,500 

4050 

Precision Landing 

384,862 

53.2 

157 kt 

400 

4069 

Lateral Offset Landing 

384,862 

53.2 

159 kt 

750 

4090 

15 kt Crosswind Landing 

384,862 

53.2 

159 kt 

1,200 

4093 

25 kt Crosswind Landing 

384,862 

53.2 

159 kt 

1,200 

4095 

35 kt Crosswind Landing 

384,862 

53.2 

159 kt 

1,200 

4129 

Category I Approach & Landing 

384,862 

53.2 

159 kt 

1,500 

5010 

Stall at Idle Power 

384,862 

53.2 

180 kt 

10,000 

5070 

Emergency Descent 

614,864 

47.2 

M 2.4 

55,400 

5080 

Low-Speed Excursion in Climb 

384,862 

53.2 

250 kt 

10,000 

5090 

High-Speed Excursion in Descent 

384,862 

53.2 

250 kt 

20,000 

7035 

One Engine Out Takeoff 

649,914 

48.1 

0 

0 

7036 

One Engine Out Cross wind Takeoff 

649,914 

48.1 

0 

0 


For some of the tasks, a formal flight card was developed to specify to the pilot what maneuver was to be performed 
within a specified tolerance of one or more parameters. These flight cards did not exist for some of the 8000 and 9000 series 
of tasks. The 8000 and 9000 series tasks, which appear in the run log listings ( Appendix E), were considered experimental 
or special-purpose and were not intended to be evaluated by a pilot, or they represented experimental variations on a task for 
which a flight card existed. 

The remainder of this appendix is composed of the formal flight cards that correspond to the tasks listed in table F-l 
above. 
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2010 Standard Acoustic Takeoff I Date: Pilot: 
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2011 Standard Acoustic Takeoff - alternate procedure I Date: Pilot: 
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19-Aug 97 EBJ 

3084 Heading Change in Supersonic Cruise - A/T off I Date: Pilot: 
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19-Aug 97 EBJ 

3086 Heading Change in Low Altitude Cruise - A/T off I Date: Pilot: 
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13 Apr 98 MDC 

3384 Heading Change in Supersonic Cruise - Small Inputs I Date: Pilot: 
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3385 Heading Change in Subsonic Cruise - Small Inputs I Date: Pilot: 
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5-Dec-96 DLR 

4090 15 Kt Crosswind Approach and Landing I Date: Pilot: 
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5-Dec-96 DLR 

4093 25 Kt Crosswind Approach and Landing I Date: Pilot: 
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5-Dec-96 DLR 

4095 35 Kt Crosswind Approach and Landing I Date: Pilot: 
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4129 IAG Straight-In Approach & Landing - Cat 1 I Date: Pilot: 
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5010 Stall at Idle Power I Date: Pilot: 
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5080 Low-Speed Excursion during Climb (subsonic) I Date: Pilot: Runs: 
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